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Abstract. By calculating the non-equilibrium parameter of the probability distribution function and the singularity spectrum 
of multifractal we have quantified the dynamical heterogeneity in strongly correlated many-body systems. 



INTRODUCTION 

We have numerically studied the dynamical heterogene- 
ity in strongly correlated many-body systems, for exam- 
ple, the critical spin statefl], the supercooled liquid near 
the glass transition[2] and the turbulence[3]. These sys- 
tems are scale-invariant and multifractal. In this paper 
we briefly review how we quantify the dynamical het- 
erogeneity and discuss in detail the turbulence simulation 
data recently obtained. 

In order to quantify the dynamical heterogeneity we 
have two measures. One is the non-equilibrium parame- 
ter q of the Renyi-Tsallis distribution function. The other 
is the singularity spectrum /(a). These two measures are 
closely related and important ingredients of the multi- 
fractal analysis. 

In the numerical study of the critical spin state[l] we 
have found that the ^-parameter represents the degree of 
non-equilibrium. We have actually quantified the devia- 
tion of the q-parameter from the equilibrium value ac- 
cording to the spatio-temporal scale of the observation. 

In the numerical study of the the supercooled liq- 
uid near the glass transition[2] we have found that the 
broken-bond distribution which reflects the nature of the 
cooperatively rearranging region is well described by the 
singularity spectrum /(oc) of multifractal. The width of 
the spectrum becomes broader as approaching the glass 
transition. Such a broadening is similar to that observed 
in the numerical study of the Anderson-localization tran- 
sition. 



TURBULENCE SIMULATION 

In the numerical study of the turbulence[3] we have 
adopted the lattice Boltzmann method. Our numerical 
turbulence on 200 3 cubic lattice points is sustained by 
a random forcing and the Reynolds number is about 500. 
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FIGURE 1. The intermittent distribution of the strength of 
the vortices. 



mean-field analysis 

Turbulence is one of the typical phenomena with mul- 
tiscale motions. Each phenomena of a scale strongly cou- 
ples with all the other scales of turbulent motion. In or- 
der to analyze such a system a scale-dependent entropy, 
the so-called e-entropy h(€), works well. For example, 
the time series of the velocity field in turbulence leads 
to a non-trivial scaling relation, h(e) °c £~ 3 , expected 
from Kolmogorov's scaling. Here e is the scale of the 
observation. The existing experimental data are consis- 
tent with this scaling. By our numerical experiment sim- 
ulating the Navier-Stokes equation we have shown the 
consistency of the scaling. In contrast to the evaluation 
of the energy spectrum, which is usually employed for 
testing Kolmogorov's scaling and determined by the two- 



point correlation function, Kolmogorov's scaling is eas- 
ily observed in the e-entropy, since it is determined by 
the mean of the exit time from the observation window 
and fluctuations are averaged out. Such a unifractal scal- 
ing is a mean-field description and fluctuations can be 
taken into account in a multifractal analysis as shown 
in the following. In the above mentioned mean-field de- 
scription the velocity difference in temporal and spatial 
directions have the same fractal scaling exponent. Thus 
we have confirmed Kolmogorov's scaling and Taylor's 
hypothesis at the same time. 

fluctuation analysis 

In order to analyze the fluctuations we use two mea- 
sures, the non-equilibrium parameter q of the Renyi- 
Tsallis distribution function and the singularity spectrum 
/(«)■ 

The probability distribution function (PDF) for the 
velocity or vorticity field is non-Gaussian. The non- 
Gaussianity is quantified by the ^-parameter. The strong 
non-Gaussianity is observed at small spatio-temporal 
scale comparable to that for the coherent vortex in the 
PDF for the velocity or vorticity difference between two 
space-time points. By using the wavelet denoising we 
have clarified that the non-Gaussianity or dynamical het- 
erogeneity results from the existence of the coherent vor- 
tex which is a strongly-correlated non-equilibrium re- 
gion. 

The strength of the vortices is shown in Fig. 1 to be 
intermittent. The spatial distribution of the coherent vor- 
tices is quantified by the singularity spectrum of multi- 
fractal /(a) as shown in Fig. 2. 
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FIGURE 2. The singularity spectrum /(a) of the spatial 
distribution of the coherent vortices. 



The remainder of this subsection is devoted to some 
speculations. As seen in our previous study [2] the width 
of the singularity spectrum depends on the degree of in- 
termittency so that we expect some Reynolds-number 
dependence of the spectrum. It has been discussed by 
many authors and should be claryfied in future system- 
atic study. In our present study the density of the vortex is 
dilute. Thus the singularity spectrum in Fig. 2 describes 
the spatial distribution of relatively free vortices. As the 
Reynolds number is increased the density increases. In 
this case the interaction among vortices becomes impor- 
tant and the correlation length of the vorticity fluctuation 
becomes large. While we can observe the vortex only as 
an individual elementary excitation in our numerical ex- 
periment, some collective excitation is expected to dom- 
inate at higher Reynolds number. In the limit of diver- 
gently large Reynolds number the correlation length be- 
comes divergently large so that we can expect full scale- 
invariance. We can find a resemblance to the case of the 
scaling theory in polymers where an ideal scaling rela- 
tion is realized for dense solutions where polymers are 
strongly entangled. In the limit of high Reynolds number 
each boxes counting the coherent vortex for calculating 
/(a) will be filled by almost equal number of vortices 
so that intermittency will disappear and unifractal Kol- 
mogorov's scaling will prevail. 

summary for turbulence simulation 

Numerical simulation data, in real space and time, for 
a forced turbulence on the basis of the lattice Boltzmann 
method have been analyzed by unifractal and multifractal 
schemes. 

Our new findings are summarized into two points. 
First in the unifractal analysis using the exit-time statis- 
tics we have verified Kolmogorov's scaling and Taylor's 
hypothesis at the same time. Second in the analysis us- 
ing the Renyi-Tsallis PDF and the wavelet denoising 
we have clarified that the coherent vortices sustain the 
power-law velocity correlation in the non-equilibrium 
state. 

Finally in the multifractal analysis it is clarified that 
the intermittent distribution of the coherent vortices in 
space-time is described as a multifractal. 
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